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Neutrinoless tau-lepton decays into either three leptons (r“ ^ or one lepton and one Ag meson(r“ —> 

£~Kg) where lepton £ means either an electron or muon, have been searched for using 48.6 fb“^ of data collected 
with the Belle detector at the KEKB e^e“ collider. No evidence for candidate decays are found in any channel. 
Therefore we set 90% confidence level upper limits on the branching fraction for 8 different decay modes. These 
limits are more stringent than those set previously and reach to the 10“^ level. 


1. Introduction 

The recent observation of neutrino oscillations 
by Super-Kamiokande[l] and SNO[2] suggests 
non-zero neutrino masses and flavor mixing in 
the lepton sector. This implies that charged lep¬ 
ton flavor violating processes such as t~ —> p,“7, 
T~ —> fjL~fifjL, T~ also occur at some 

level. It is therefore important to search for lep¬ 
ton flavor violation in rare decays of charged lep¬ 
tons. 

In the standard model(SM), lepton flavor vi- 
olating(LFV) decays of the charged lepton are 
highly suppressed even if we consider the effect 
of neutrino mixing[3]. However, many extensions 
of the Standard Model predict LFV decays of 
changed leptons [4]. For example, in the min¬ 
imum SUSY standard model with right-handed 
neutrinos[5,6] or in a SUSY grand unified theory 
model[7], flavor mixing in the slepton mass ma¬ 
trix at the TeV scale can induce a LFV transition 
T —> /i(e). The predicted rate of r —s- ^7 decays 
in these models depends on the model parameters 
but ranges from 0(10“®) to 0(10“®) , which is in¬ 
terested area of the current and future H-factory 
experiments. On the other hand, other models, 
with multi-Higgs [8], extra Z gauge bosons [9], 
R-parity violating interaction[10] or generic mod¬ 
els involving heavy neutrinos [11,12], predict rela¬ 
tively large branching fractions for r“ ^ 


decays. The maximum allowed value can reach 
the 0(10“^) - 0(10“®) level. Searches for flavor 
violating decays in various decay modes are thus 
important and can provide strong constraints on 
possible new physics. 

Our recent results for the decay t“ ^ /i“7 are 
reported in another paper[13]. In this paper, we 
present the results of a search for LFV decays into 
three leptons or one lepton and one K'g meson; 

T“ ^ (1) 

r“ ^ i-Kl (2) 

where £i stands for either an electron or muon[14]. 

In previous experiments, the most stringent up¬ 
per limit on the branching fraction for the decay 
r“ ^ has been set by the CLEO exper¬ 

iment. They obtained an upper limits of (1-2) 
X 10“® at 90% confidence level from a data sam¬ 
ple of 4.79 fb“^[15]. For the decays t“ ^ (-~Kg, 
previously published upper limits come from the 
Mark H Collaboration and are of 0(10“^)[16]. 
The CLEO group has recently reported the im¬ 
proved upper limits iJ(r“ ^ e~Kg) < 9.1x10“^, 
B{t~ ^~Kg) < 9.5 X 10“^ at 90% C.L. using 
a 13.9 fb“^ data sample[17]. 

2. Data and Event Selection 

The data used in this analysis were collected 
with the Belle detector[18] at the KEKB[19] 
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asymmetric energy e“'"e“ collider operating at a 
center-of-mass energy ^/s ~ 10.6 GeV. Among 
the recorded samples, we use 48.6 fb“^ of data, 
corresponding to 44.2 million r-pair productions, 
taken in the period between 1999 and 2001. 

The Belle detector is a general purpose detec¬ 
tor having excellent capabilities for precise vertex 
determination as well as the particle identifica¬ 
tion. Tracking of charged particles is performed 
by a three layer double-sided silicon vertex detec¬ 
tor (SVD) and a 50 layer cylindrical drift chamber 
(CDC) within the 1.5 T magnetic field. Charged 
hadron are identified by means of dE/dx from 
the CDC, signal pulse-heights from the aerogel 
Cerenkov counters (ACC), and timing informa¬ 
tion from the time-of-flight scintillation counters 
(TOF). Muons are detected by a 14 layers of re¬ 
sistive plate counter interleaved with iron plates 
(KLM). Photons and electrons are detected using 
a CsI(Tl) electromagnetic calorimeter (ECL). 

In order to determine the selection criteria, a 
Monte Carlo (MC) studies are performed using 
the KORALB/TAUOLA program [20] for the r- 
pair generation, the QQ generator[21] for BB 
and qq continuum processes, the AAFHB[22] pro¬ 
gram for two-photon processes, the BHLUMI[23] 
program for radiative Bhabha events and the 
KKMC[24] program for radiative /r-pair events. 

Regarding the signal r-pair sample, one r de¬ 
cays into or assuming 3- or 2-body 

phase-space distribution, respectively, in the rest 
frame of r, and the other r decays generi- 
cally. The detector response is simulated by a 
GEANT[25] based program. The MC events thus 
prepared are passed through the same reconstruc¬ 
tion and analysis program as for the actual data. 
Based on the study of simulated signal and back¬ 
ground samples, the candidate signal events are 
required to satisfy the criteria described below. 

2.1. Topological Event Selection 

To attain a large detection efficiency, we re¬ 
quire 1-prong decay {B{t “1-prong”) = 85.4%) 
for one t lepton in the reaction e“'"e“ ^ 
and the other t decays into three charged final 
state. Accordingly the signal candidate events 


must contain four charged tracks with zero net 
charge. The tracks should be well-reconstructed 
and have transverse momentum Pt > 0.1 GeV/c 
within a polar angle 17° < 6iab < 150° in the 
laboratory frame. The tracks include the ones 
originating from the beam interaction regions as 
well as the daughter tracks from the secondary 
vertex. The latter tracks are added in order 
to include the decay products of the Kg meson. 
To reject beam-gas events, the distance of closest 
approach of each non-ATg track to the interaction 
point must be within ±1 cm transversely and ±5 
cm along the beam. Photons are defined as neu¬ 
tral ECL clusters with an energy > 0.1 GeV. 
Photons must be separated by at least 30 cm from 
projection of any charged tracks on the surface of 
the calorimeter. 

We divide each event into two hemispheres in 
the e“'"e“ CM system with the plane perpendic¬ 
ular to the thrust axis. The thrust axis is cal¬ 
culated from the momenta of charged tracks and 
photons. In each hemisphere, one(tag) side must 
contain one charged track. The other(signal) 
hemisphere must contain three identified leptons 
or one identified lepton and one Kg meson. There 
is no explicit cut on the maximum number of pho¬ 
tons on both hemispheres in order to maintain 
high detection efficiency while minimizing the de¬ 
pendence on MC simulation of showers. 

In order to reject the two-photon processes, 
total transverse momentum in the event is re¬ 
quired to be larger than 0.5 and 0.2 GeV/c for 
T ^ and r ^ £Kg decays, respectively. 

To reduce the radiative Bhabha and radiative /t- 
pair background, the collinearity angle 6coi be¬ 
tween signal and 1-prong side is required to sat¬ 
isfy 6coi < 175°. 

2.2. Lepton Identification 

Lepton identification is substantial in this anal¬ 
ysis. Electrons are identified by means of a likeli¬ 
hood ratio combining information on the ratio of 
the cluster energy in the ECL to the track mo¬ 
mentum measured in the CDC and dE/dx in the 
CDC. The electron identification efficiency is es¬ 
timated to be 92% for the momentum range be- 
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tween 1.0 GeV/c and 3.0 GeV/c in the labora¬ 
tory frame. The rate for misidentifying hadrons 
as electrons is quite small (0.3%). In order to 
correct for the energy loss from the radiation in 
the material of the SVD region, the energy for 
the electron candidate is recalculated by adding 
the energy of radiated photon clusters, if an EGL 
cluster with energy less than 1.0 GeV is detected 
within a cone angle of 10° around the flight di¬ 
rection of the electron candidate track. 

Muon are required to have a well reconstructed 
track in the muon system, comprised of 14 lay¬ 
ers of iron plates interleaved with the KLM. The 
muon probability is calculated from two variables: 
the difference between the range calculated by the 
momentum of a particle and the range measured 
by the KLM as well as the of the KLM hits 
with respect to the extrapolated track. The iden¬ 
tification efficiency is estimated to be ~91% for 
momenta in the laboratory frame greater than 
1.2 GeV/c. The rate for misidentifying hadrons 
as muons is 2.0%. 

2.3. Kg Selection 

Kg candidate is selected from candidates 
assuming their decay products are charged pi¬ 
ous. The candidates are reconstructed from 
all pairs of oppositely charged tracks with the fol¬ 
lowing requirements; 

i) the closest distance between two helices of the 
candidate pair must be less than 1 cm in the 
coordinate along the beam axis, 

ii) the flight length from the beam interaction 
point to the crossing point in the transverse 
plane must exceed 2 mm, 

iii) the angle between the direction of crossing 
point and the direction of the momentum 
must be less than 10°. 

Figure 1 shows the tt+tt” invariant mass dis¬ 
tribution for the data after applying topological 
conditions for events which include four charged 
tracks. A clean Kg signal is seen with very low 
level background. The solid curve in the figure is 
a fit with a double Gaussian for the signal and a 
linear function for background. The fitted mean 
of the Gaussian fit is consistent with the nominal 


Kg mass. The mass resolution for the Kg signal 
is 5.1 ± 0.1 MeV. The Kg candidates must have 
an invariant mass within three stan¬ 

dard deviations of the nominal Kg mass, 0.485 
< M^+^- < 0.510 GeV/c2. 

2.4. Signal Selection 

To identify the neutrinoless tau decays, a signal 
region is defined by using two kinematical vari¬ 
ables; 

A ZT** — Z?* Z?* 

l\rj — ^signal -^beam 

AM = M,,gnal-Mr, 

where and Mr are the beam energy in the 

e+e“ CM system and the r lepton mass, respec¬ 
tively. and Msignai are the sum of the 

energy in the signal side and its invariant mass, 
respectively. is measured in the e+e“ CM 

system. Figure 2 shows the AE* vs. AM plots 
for the signal Monte Carlo sample. The signals 
are concentrated in the vicinity of AE* = 0 and 
AM = 0. The tail of AE* on the lower side is 
caused by initial photon radiation in the r-pair 
production process. A small tail, seen on lower 
side of AM for the processes involving electrons, 
is caused by the leakage of electron showers in 
EGL. The signal region for the LFV decay modes 

-0.29 < AE* < 0.11 GeV 
-25 < AM < 25 MeV/c^ 

is shown by the dashed lines in Fig. 2. These 
requirements correspond to three standard devi¬ 
ation limits. Figure 3 shows the AE* vs. AM 
plots for the experimental data after applying all 
selection cuts except for AE* and AM. Each de¬ 
cay mode includes the decay products of the t~ 
lepton as well as the charge conjugate decay mode 
from the lepton. 

2.5. Background and Efficiency 

The backgrounds remaining in the sample can 
be estimated from the side-band region shown 
in Fig.3. The background is negligible for the 
modes with r decays into three leptons, while 
somewhat larger background is observed for the 
decays t~ ^~Eg. Study of the background 
MG indicates that this background comes mainly 
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from the decays t~ Tr~ {K~)Kgi'r, where a 
charged pion or kaon is misidentified as an elec¬ 
tron or muon. 

In order to estimate the background in the sig¬ 
nal box, a fit with a linear function is carried out 
for the AE* distribution of the side-band region 
for the sample -25 < AM < 25 MeV/c^. Since 
the statistics is low, the slope of the function in 
the fit is fixed to the one determined by the fit 
of the same distribution but selected with looser 
or no lepton identification requirements *. The 
results for the estimated backgrounds are shown 
in Table 1. The error on the background comes 
from uncertainty in the fitting procedure. 

The signal detection efficiency is determined 
from the signal MC. The efficiencies for each of 
the LFV decay modes are listed in Table 1. The 
efficiency is about 10% for t~ £1(2^3 de¬ 
cays, and is about 16% for t~ decays. 

The latter value includes the branching fraction 
for Kg 7r+7r“ (68.60 ± 0.27%) as well as the 
branching fraction of the 1-prong decay of the tau 
lepton. 

The systematic errors for the efficiency deter¬ 
mination come from tracking (2%/track), Kg se¬ 
lection (0.4%), electron identification (1.5%) and 
muon identification (1.3%). The uncertainty in 
the luminosity measurement is estimated to be 

1.5%. In the calculation of the upper limit on the 
branching fraction, these systematics are taken 
into account by reducing the efficiencies by Ict 
from its normal MC values. 

3. Results 

No candidate decays are observed in the signal 
box for any of the eight decay modes, as can be 
seen in Fig. 3 and Table 1. 

Using the Feldman and Cousins method [26, 
27] for zero observed events and the estimated 
number of background events in the signal box, 
upper limits are calculated for each decay mode. 
The final results are summarized in Table 1 for 
the 8 decay modes. 

* Since the background are mainly from the fake lepton 
identification, its shape is expected to be similar in all 
cases. 


The upper limits on the branching fractions are 
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at 90% confidence level. 

The limits obtained in this analysis for the de¬ 
cays T~ £ii2^3 are approximately one order 
of magnitude more restrictive than the current 
limits. The upper limits for t~ —> £~Kg decays 
are a factor of 4-5 more stringent than recently 
reported by the CLEO experiment [17]. 

Since the background level is still low, we ex¬ 
pect to be able to search for the LFV decay mode 
at down to the 0(10“®) level in the near future 
at high luminosity B-factory experiments. 
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Table 1 

Summary of detection efficiency, observed events, expected background and 90% C.L. upper limits on 
the branching fractions. 


Mode 


EfRciency(%) 

Observed 

events 

Expected 

background 

Upper limit 
(This Exp.) 

Upper limit 
(PDG2002) 

r~ 

e e+e 

12.1 ±0.7 

0 

0.4 ± 0.3 

2.7 X 10"' 

2.9 X lO"*’ 

T~ 


9.0 ±0.6 

0 

0.0 ± 0.1 

3.8 X 10-" 

1.9 X 10"*’ 

T~ 


10.4 ±0.6 

0 

0.5 ± 0.4 

3.1 X 10"^ 

1.8 X 10"*’ 

T~ 

/i“e+e“ 

11.9 ±0.7 

0 

0.8 ± 0.4 

2.4 X 10"^ 

1.7 X 10"^ 

T~ 

li~ li~ 

10.8 ±0.6 

0 

0.0 ± 0.0 

3.2 X 10"^ 

1.5 X 10"^ 

T~ 


12.2 ±0.7 

0 

0.0 ± 0.0 

2.8 X 10"^ 

1.5 X 10"^ 

T~ 

e-K'^ 

16.0 ±0.7 

0 

1.7 ± 0.5 

2.9 X 10"' 

1.3 X 10"^ 

T~ 

IFKf— 

17.4 ±0.7 

0 

2.0 ± 0.5 

2.7 X 10"^ 

1.0 X 10"^ 



Figure 1. 7r“''7r invariant mass distribution for the experimental data sample. The solid line is the result 
of a fit with a Gaussian for the signal and a linear function for the background. 
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Figure 2. AE* vs AM plots for the LFV decay modes in the signal Monte Carlo sample. The dashed 
lines indicate the signal region. 
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Figure 3. AE* vs AM plots for the experimental data sample. The charge conjugate decay mode is also 
includes. The dashed lines indicate the signal region. 





















































